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The transport properties of quantum dot (QD) systems based on double-walled carbon nanotube 
(DWCNT) are investigated. The interplay between microscopic structure and strong Coulomb inter- 
. . . action is treated within a bosonization framework. The linear and nonlinear G-V-Vg characteristics 

' of the QD system is calculated by starting from the Liouville equation for the reduced density ma- 

\ trix. Depending on the intershell couplings, an 8-electron periodicity of the Coulomb blockade peak 

■ spacing in the case of commensurate DWCNT QDs and a 4-electron periodicity in the incommensu- 

(N . rate case are predicted. The contribution of excited states of DWCNTs to the nonlinear transport 

is investigated as well. 
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I. INTRODUCTION 



After being discovered in lOOl'^, carbon nanotubes (CNTs) have been widely used in nano-devices because of their 
unique properties^i^. CNTs may be either single-walled (SWCNT) or multi-walled (MWCNT) depending on the 
' number of graphene sheets wrapped into concentric cylinders. Due to the quasi-one-dimensional characters of their 
CZ5 \ electronic structures, long SWCNTs exhibit Luttinger-liquid behavior— i^i^ii^. SWCNT quantum dot (QD) systems 
1^ • have also been fabricated, which consist of finite length SWCNTs weakly connected to the source and drain leads 
^ ■ and capacitively coupled to a gate electrod o^^d^d^d^ ^ p^^ ^aw bias, the SWCNT QD systems show Coulomb blockade 
^ . behavior because of the strong Coulomb interactions in the QDs and the poor transparencies of the contacts between 
' the QD and the leadsi^. Because of the short lengths of SWCNTs, the addition energy needed to add an extra 
^ , electron to the QD depends on both the Coulomb interaction and on the energy level spacing. Unlike the traditional 
I ' two-dimensional semiconductor QD systems with irregular Coulomb blockade patterns, which have to be understood 
^ \ statistically^, QDs based on SWCNTs show regular Coulomb blockade patterns, which originate from the regular 
Q • electronic structure of the SWCNTs. Because of the spin degeneracy of two bands crossing at the Fermi points in 
Q \ metallic SWCNTs, the stability diagrams of SWCNT QD systems exhibit a 4-electron periodicity of the Coulomb 
'— diamond size a"'^°d^d^d3,i6 _ ^pj^g stability diagrams of the SWCNT QD systems have been explained by using the 
^ ■ mean- field theory developed in Ref. [ij which includes a noiizero exchange energj*i2iiiii^iH. Recently, the energy 
spectrum of SWCNT QD has been calculated in Refs. 18' and 1J| beyond mean field. For QD systems with moderate- 
\^ ' to-large radius SWCNTs, the exchange energy can be ignored^-, and the stability diagrams can also be quantitatively 
, explained within a bosonization approac h^^d^ . By suitable choice of parameters these theories can reproduce the 
■ same low bias spectra of SWCNT QDs, and only the excitations measured at high bias are predicted differently by 
' a mean-field approach or by a bosonization method because of the different treatment of the Coulomb interactioi*i^. 
. Although the excitations of SWCNT QDs have already been measured^^, the quality and the range of the measured 
' excitations cannot be used to determine the validity of these two methods and further experiments are needed, 
r ■ ■ So far, the properties of MWCNT QD systems have not been fully explored^ The experiment in Ref. '2T showed 
' that the stability diagrams of MWCNT QD systems have a 4-electron periodicity of the Coulomb diamond sizes. 
J> ', The simplest MWCNT QD is the one formed by a double- walled carbon nanotube (DWCNT), which consists of two 
r ■ concentric shells. Depending on the ratio between the unit cell lengths of the two shells, a DWCNT may be either 
rN ' commensurate (c- DWCNT), if the ratio is a commensurate number, or incommensurate (i-DWCNT), if the ratio is an 
incommensurate number. It has been shown that the effective intershell coupling depends on the type of DWCNTs. 
■ " " '■ At low energies, that is, near the Fermi energy, the effective intershell coupling is negligible in i-DWCNTs but large in 
c-DWCNTs while it cannot be ignored in both type of DWCNTs at high energies2Si2^i^i2^. Both types of DWCNTs 
with long lengths can be described by Luttinger liquid theory when Coulomb interactions are included^^'^^. Because 
of their intermediate-to-large radii, we expect that the exchange energy may be ignored in DWCNTs. Therefore, the 
bosonization approach, which includes forward scattering processes exactly, can be used to describe the properties of 
DWCNT QD systems as well. 

In this paper, we consider a QD system formed by a finite length DWCNT with two metallic shells, where we 
include all forward scattering processes. The bosonization approach enables exact diagonalization of the interacting 
DWCNT Hamiltonian. Finally, the linear and nonlinear transport properties of the system are investigated by solving 
the Liouville equation for the reduced density matrix to lowest order in the coupling to the leads. 

The paper is organized as follows. In Sec. HIl the Hamiltonian of a DWCNT QD system is derived. The energy 
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FIG. 1: (Color online) Schematic experimental setup of a double- walled carbon nanotube (DWCNT) quantum dot (QD) system. 
A finite length DWCNT is deposited on a substrate and weakly connected to the source and drain leads through its outer shell. 
A gate electrode is capacitively coupled to the DWCNT QD and controls the electrochemical potential in the QD. The dashed 
lines denote the inner shell in the DWCNT. 



spectrum of a finite length DWCNT with strong Coulomb interactions and open boundary conditions is then obtained. 
Transport properties of DWCNT QDs are calculated in Scc. lIIII The results for the linear and nonlinear conductances 
of both c- DWCNT and i-DWCNT QDs are presented in Sec. |TVl Finally, the conclusion is drawn in Sec. |Vl 



II. MODEL AND METHOD 



As schematically shown in Fig. [U the QD system consists of a DWCNT with two metallic shells deposited on a 
substrate. The source and drain leads are connected to the outer shell of the DWCNT. The segment of the DWCNT 
(of about several hundred nanometers long) between two leads forms a QD. A gate electrode is capacitively coupled 
to the QD and controls the electrochemical potential in it. As we are only interested in the Coulomb blockade regime, 
we assume that the QD is weakly contacted to two leads, that is, the transparencies of the contacts are very poor 
and the conductance of the QD system is much smaller than the conductance quantum 2e^//i. The Hamiltonian of 
the whole system can be separated into several parts, 

H = -fflcads + Hqd + Ht + Hg, (1) 

where Hqd is the Hamiltonian of the QD system and its explicit form will be derived in the following subsection. The 
source (s) and drain (c?) leads are described by Fermi gases of non-interacting quasi-particles and the Hamiltonian of 
the leads is 

i^ieads = 51 (^'k - ey/)c|k„qk^, (2) 

l—s,d kcr 

where e is the elementary charge and Vi is the voltage in the lead I. The operators cj^^^ and Ci^^ arc the creation and 
annihilation operators of a quasi-particle with wave vector k and spin cr = ± in the lead I. The Hamiltonian of the 
gate is 

Hg = -e^gAf, 

where /Xg is the chemical potential in the gate and the operator J\f accounts for the total number in the QD system. 
Ht is the tunneling Hamiltonian describing the tunneling between the QD and the two leads and it has the form 

Hr^Y^Hj '^^ Twiy)^lM<^ia{v) + H.c, (3) 

l—s.d /3cr 

where ^laiy) = X^k '^k(r)C[j.;^ is the electron annihilation operator in the lead / and ^'^^^.(r) is the electron operator 
in the shell j3 whose explicit form will be given in Sec. IH Al 



A. Low energy non-interacting Hamiltonian of DWCNT 



In general, the energy spectrum of a SWCNT or of a DWCNT without electron-electron interactions can be obtained 
by using a tight-binding model for the pz orbitals in carbon atoms^. In particular, we shall view in the following a 
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DWCNT as two tunneling coupled SWCNT shells. We denote with the index /? = ± the outer/inner SWCNT shell. 
A metallic SWCNT shell (3 within periodic boundary conditions has two independent Fermi points {^K^ p). Their 
positions depend on the chirality of the shell and in general are different for different SWCNT shells. At the Fermi 
points, the lowest conduction and the highest valence bands touch each other as shown in Fig. [SJa). As the next 
conduction and valence bands are separated by a large gap (about 1 eV)^, we will only consider the lowest conduction 
band and the highest valence band in our calculations. The energy dispersion near the Fermi points is linear—, see 
Fig. IHa), and is given by 

£R/L{li) = ±flVFK, (4) 

where the wave vector k is measured with respect to the Fermi points and the Fermi velocity in SWCNTs is vf ^ 
8 X 10^ m/s. Hence, at each Fermi point, there are two branches corresponding to the right (+) and left (— ) moving 
electrons. The Bloch waves for the electrons in these branches in a shell P — ± are 

^PprFnir) = e"-'''(pf}rF{r), (5) 

where r = {u,v) and u and v are along the nanotube axis and the circumference directions, respectively (cf. Fig. [3]). 
The index r = R/ L = ± denotes the right and left- moving electrons. The periodic function is 

V>/3rF (r) = — i= V e'^-^fpprF x(r - R - Tp) , (6) 

where the index R denotes the lattice vector of the graphene sheet, is the number of carbon atoms in the shell 
f3, p is the index for the two graphene sublattices and Tp is the vector giving the positions of the two different atoms 
in a unit cell. The index F is for the Fermi points in the shell and F denotes the Fermi points in a graphene sheet. 
They are related as F = F • u with u the unit vector along the nanotube axis. The coefffcients /'s depend on the 
chirality of the shell (m, n) as^ 

fpAvF = (-^sgn(Fr)(ri + m) + ^sgn(r)(m - n)) , (7) 

fpBrF = (^) 



where C = + mn + and the function x(r — R) is the Pz orbital wave function. Because we consider a finite 
length shell, we have to use the open boundary condition (OBC) instead of the periodic boundary condition along 
the tube axis (cf. Fig. [5]). The wave function in the shell f3 satisfying the OBC has the forroi^ 

= -^{fPR/LKoni^) " ^0L/R.-Ko-^i^))> (9) 

and the wave vectors k are quantized as 

K=^{m, + Ap), m« = 0,±1,±2,--- , (10) 

where L is the length of the nanotube. The mismatch of the Fermi points is < Ap — Ko^pL/n — [Xo./jL/tt] < 1, 
where [• • • ] gives the integer part of its argument. The Hamiltonian of a finite length non- interacting shell (3 is thus 

^° ^^^&^(^)f^^FKclf^^Cf^f^^, (11) 

where r = R/L — ± is the index for the left and right moving electrons with the OBC. The operators cj^^^^ and Cfj^^^ 
are the creation and annihilation operators of an electron in the branch f with the wave vector k and spin a. 

Let us now see how the spectrum gets modified when looking at DWCNT. At low energies, the intershell couplings 
are different for c-DWCNTs and i-DWCNTs. The intershell couplings in i-DWCNTs are negligible while they are quite 
strong in c-DWCNT o^^'^^i^^'^'^i^^i^^ . Therefore, the non-interacting Hamiltonian of an i-DWCNT is the combination 
of the Hamiltonians of two shells, 

-ff^DWCNT = sgn(r)nwi.Kc)3^^^c^f;^^ 

3r(TK 

(12) 

= sgn(f)(m«eo-f^ A0eo)4^«^c^f.^<„ 
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FIG. 2: (Color online) Energy spectra of metallic single-walled (SWCNT) and double-walled (DWCNT) carbon nanotubes 
with periodic boundary conditions (PBC). (a) Energy spectrum of a metallic SWCNT. There are two Fermi points and two 
branches L/R with left/right- moving electrons at each the Fermi point, (b) Energy spectrum of an incommensurate DWCNT 
(i- DWCNT). It consists of the energy spectra of the outer and inner graphene shells (±), which are not coupled to each other 
because of the vanishing intershell coupling, cf. Eq. (|12|) . (c) Energy spectrum of a commensurate DWCNT (c-DWCNT). 
Because of the finite intershell coupling, it is composed of the bonding and anti-bonding bands (±), which are shift vertically 
along the e-axis, cf. Eq. (fT5)l . 



where eo — HvpTr/L is the level spacing and we have used the quantization relation for k, Eq. (fTO|) (cf. Fig. [4l^b)). 
On the other hand, the non-interacting Hamiltonian of a c-DWCNT contains also the contribution from the intershell 
coupling^, 



c-DWCNT = XI XI Sgn(^)?lWFK4-^^C^^„^ + X X ^4r«<TC/3'fK^ + H-C, 



(13) 



(3 rcTK 



0(3' rcTK 



where t is the intershell coupling and we assume that it is a constant in the low energy regime. The Hamiltonian 
Eq. (jl3p can be diagonalized by using the bonding and anti-bonding basis. 



(14) 



where v — ± is the index for bonding and anti-bonding states, respectively. The non-interacting Hamiltonian of a 



FIG. 3: (Color online) Cross section of a DWCNT. Atoms A and B in two shells of radii R+ and R-, respectively, are projected 
onto this cross section. Such atoms are described by the coordinates {u+,v+) and {u-,V-), where u± are along the tube a^cis 
and v± measure the atom positions on the outer/inner circumference. 






FIG. 4: (Color online) Energy spectra of metallic single-walled (SWCNT) and double-walled (DWCNT) carbon nanotubes with 
open boundary conditions (OBC). (a) Energy spectrum of a metallic SWCNT. There are two branches L/R with left/right- 
moving electrons. The parameter eo is the level spacing and A describes the mismatch of the Fermi point, cf. Eq. (|10[l . 
(b) Energy spectrum of a commensurate DWCNT (c-DWCNT) and (c) of an incommensurate DWCNT (i-DWCNT). The 
parameter A± describes the mismatch of the Fermi point in the shell ± and for a c-DWCNT is A+ — A_ = A. The parameter 
describes the mismatch of the states in two bands, cf. Eq. (I16p . 



c-DWCNT in the new basis becomes 

-"c-DWCNT - 



X! (sgn(^)("^Keo + Aeo) + sgn(!/)C£o)4 



(15) 



where for c-DWCNT is = A and the parameter C is defined as 

< C - Veo - [t/eo] < 1, 
which describes the mismatch of states in two bands (cf. Fig. IHc)). 



(16) 



Coulomb interaction Hamiltonian of DWCNTs 



In quasi-one dimensional electronic structures as CNTs, Coulomb interactions are not fully screened and can strongly 
influence the properties of CNTs^i^i^i^i^ i^^i^^ . The total Coulomb interactions in an i-DWCNT can be expressed by 
the following Hamiltonian, 



Tjint 

-"i-DWCNT 



\ jj dvidv2^\,{ri)-^l,^,{T2)Upii,{vx^T2)-^p,,,{v2)^pAri), 

13(3' aa' 



(17) 
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where = {ui,Vi) and u and v are along the tube axis and the circumference direction, respectively (cf. Fig. [3]) . 
The intrashell interaction is given as 

Uf30{n - V2) = , (18) 



and the intershell interaction is 



U+^{v^ - r2) = =, (19) 

(lii - U2)2 + 4i?+i?_ sin2(ui/2i?+ - V2I2R-) + Ai?2 

where e is the dielectric constant, az is the "thickness" of a graphene sheet and the distance between two shells is 
Ai? = |i?+ — R-\. The electron operators for i-DWCNTs are defined as 

^p,{v) = Y,v'0f,^{^)cpf., (20) 

fq 

Using Eq. (O we can define the ID electron operators describing the slowly varying part of the electron operators 
^',3^(r) as 

i:p,F^{u) = 4f E e'^*^"(^)''"c^,.,, (21) 



2^ , 

in terms of which the electron operators can be written as 

«'/3a(r) = ^(/J^i^^(r)c^,,^ 

(22) 

= VIE^Sn(-F)'/'/3sgn(F)fF(l")'0/3fFa(w)- 
fF 

For SWCNT shells with diameter larger than ~ 1 nm, we can only keep forward-scattering (or density-density) 
processes22i, such that the interacting Hamiltonian in an i-DWCNT becomes 

^^i-DWCNT = ^EEEE // duidu2 PpfFa{ui)Vp^f {ui - U2) P0'f' F' a' {U2) , (23) 
1313' fr' FF' aa' 

where Pi3fFaiu) = '0^fFcr('^)^/3fFCT('") electron density operator and V^^^/ is the effective one-dimensional 

Coulomb interactions given by 

V;^^/ {m,U2) ^ J2 [[ dv^dv2\x{ri-Ri)\'Uf3p'ir,~r2)\xir2~ 112)1'. (24) 

R1R2 

Let us now turn to c-DWCNTs. The total Coulomb interaction Hamiltonian of a c-DWCNT has a similar form as 
Eq. ()17|) . However, we have to rewrite it in the basis of the bonding/anti-bonding states using the transformation 
Eq. (fTi]) . In this basis, the Coulomb interaction of a c-DWCNT is given by 

K%WC^T^l E jj dv,dV2^lAv{)^l^,{v2)ij,.-{v,-V2)^,.A^2)^.M, (25) 

where = ± is the index for bonding and anti-bonding bands and the new interactions are 

= = 2U^_ = ^{U++ + C/__ + C/+_), (26) 

where [/g/3' is defined by Eqs. and The electron operators in c-DWCNTs are defined as 

rq 

= ^^^S^{F)(p^^^n(F)fFij:)4'vf<yF{u), 



(27) 
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where (p'^^q^ is a linear combination of f'^^J" , namely, (^^^*-^(r) = (^{ip^?^ (r) + sgn{h') ip'^f^ {r)^ /a/2. By using these 
electron operators and keeping only the relevant forward-scattering processes, the Coulomb interaction Hamiltonian 
of a c-DWCNT becomes 

■f^c-DWCNT = ^ X! X] // duidu2Pt,rFcriui)V^l/iui-U2)p^'r'F'a'iu2), (28) 

uv' rr' FF' go' 

where PufFaiu) — ^lfp^{u)^l>^-p^{u) is the density operator and V^^/ is the effective one-dimensional Coulomb 
interaction, 

K!/{ui,U2)^-^ [[ dv,dv2\x{ri- ■R,)\^U.Ari~r2)\xir2 -112)1' ■ (29) 

rLirL2 



C. Bosonization 



The low energy Hamiltonian of a DWCNT is the combination of the non-interacting and interacting Hamiltonians 
and it can be diagonalized by the bosonization metho d^^'^^i'^^'^'^ . First, we introduce the bosonic operatorsiSi^i^ 

, _ {pcfqal^, for i-DWCNTs, 

forc-DWCNTs, ^''^^ 

where q = -KUq/L with rig an integer. The index a = ± denotes the bonding/anti-bonding states in c-DWCNTs and 
outer/inner shells in i-DWCNTs and we will keep this convention in the rest of the paper. The bosonic operators 
obey the bosonic commutation relation 

By using these bosonic operators, the Hamiltonian of a DWCNT QD can be separated into its fermionic and bosonic 
parts, Hq£) = Hf + H^. The fermionic part Hf describes the ground state and the fermionic excitations in the 
DWCNT. The fermionic Hamiltonian of a c-DWCNT is 

^f/x-DWCNT = \^^^lra + A£osgn(f )7V;f ^ -l- (aC^^o - \^e^Mo,fa + Hf\ (31) 

OLVa 

while the i-DWCNT fermionic Hamiltonian has the form 

i?/a-DWCNT = Y \^^^lf-y + Aa£osgn(f)7V;f^ - ^eoA/'afo- + llf\ (32) 

<y.r(7 

where B.^^ is due to Coulomb interaction having the form 

olol' rcr r'cr' 

with the interaction strengths 

WSr' =^jj duldU2 V:il{u, - U2). 

Therefore, the fermionic Hamiltonian of a DWCNT QD is described by the constant-interaction model^^. 

The bosonic excitations of a DWCNT QD are described by the Hamiltonian Hf,, which can be expressed in terms 
of the bosonic operators as 

1 ca' t t ^^^^ 

q>0 cxa' rf' acr' 



8 



with the interaction strengths 



^raa 
qq 



1 
l2 



duidu2 V^[f {ui - U2) cos(gui) cos{qu2)- 



(35) 



In order to diagonalize the Hamiltonian i/f,, we need to introduce new bosonic operators ftj^^g's, where j = c, s denote 
charge/spin modes and the remaining indices S = ± and ^ = ± define total/relative modes with respect to the branch 
and shell (or bonding/anti-bonding state) degrees of freedoms, respectively. The new bosonic operators are related 
to the bosonic operators b 



i(f)g(T 



s ai 



.25.34 



(36) 



where the coefficients are given by 
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(37) 



with 



smf 



w- 



W, 



++ 



w:. 



WqY - Wq-q- 



1/2 



The other coefficients are 



Sjs^q = 1 and Cjs^q = 



(38) 



in the cases (jS^) = (c— ±), (s±±). For the total and relative charge modes (c+±), the two coefficients are interaction 
dependent 



Sc+±q 



£0 



£c+±{q) 



l £c+±{q) 
£0 



Cc+±q 



£0 



l£c+±{q) 



£c+±iq) 



where the energies of the total and relative charge modes are 



£c+±{q) =£q\/1 + SW^q'^/eo 



(39) 



(40) 



The interactions do not affect the 6 "neutral" modes, {jS£,) — (c — ±), (s ± ±) and their energy dispersions are the 
same as for the non-interacting system. 



By using the new bosonic operators, the excitation Hamiltonian can be diagonalized to be 



and the eigenstates are 



q>0 jSi 



N,m).nnT^(«kj'""'iN,o) 



q>0 jS^ 



(41) 
(42) 

(43) 



where N = {Nafa} defines the number of electrons in each of the eight branches (ara) and m — {mjg^g} describes 
the configuration of the bosonic excitations in each of the eight modes (jS^). The state |N, 0) contains no bosonic 
excitations and describes the ground state or the fermionic excited states. 
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III. DYNAMICS OF THE QD SYSTEM 

The transport properties of the DWCNT QD system can be obtained by investigating the dynamics of its density 
matrijs^. In this section, we briefly show how to derive the equation of motion for the reduced density matrix of the 
DWCNT QD system. By solving these equations we obtain the stationary current through the DWCNT QD system 
when a bias voltage is applied. 



A. Equation of motion for the reduced density matrix 

As we consider a very weak coupling between the DWCNT QD and the two leads, the tunneling Hamiltonian can 
be treated as a perturbation and we can obtain the equation for motion for the density matrix in the interaction 
picture as^ 

^^^-^-[HktlpUt)l (44) 

where pl^i {t) is the density matrix of the whole system (including the DWCNT QD and the leads) and the tunneling 
Hamiltonian in the interaction picture is 

This equation can be solved formally as 

pLit) = pLih) -Jildh [HUti),pUti)]- (46) 
Substituting the above expression of pj^^ (t) back to Eq. ([i^ , we have 



it) - ^l[Hi^{t),pUto)] + {jS ly^' t^^^*^' [^T(ii),pL(^i)]]- (47) 



dt 

As we are only interested in the transport through the DWCNT QD, we will focus on the reduced density matrix of 
the QD which is obtained by tracing out the degrees of freedom of the leads, 

/ = Tr,eads{pLt}- (48) 

Because the leads are very large comparing with the QD and the tunneling events between leads and the QD are rare, 
the effect of the QD on the leads can be ignored and the leads can be described as reservoirs remaining in thermal 
equilibrium. We use the ansatzi^, to factorize the total density matrix (t) , 

pLM - pLdsp'it) = piPdp\t). (49) 

where the density matrix of the leads, picads, is time independent and is described by the thermal equilibrium distri- 
bution, 

J. g-|3iHs/d-^l■s/d^fs/d) 

P'^/d' " rp^|g-/3(lf,/d-Ms/dA/',/d)}' 

where p,s/d is the chemical potential of the source/drain lead and (3 = l/fc^T. We further simplify Eq. (I47p by 
introducing the Markov approximation, that is, we assume that plgi{t) only locally depends on pfoj(i) and we can 
replace pL(^') by pL(0- 

We make the further assumptions that the elements of the reduced density matrix between two states with different 
charges vanish, and that the elements between two non-degenerate states with same charges also vanis h^^i^^ . Finally, 
the master equations of the reduced density matrix can be expressed in Bloch-Redfield form^^i^ 

p^'^"(o--E^n™.."^l.^"W+ E EE^f^fc^'^/^i.^'-'W' (50) 

kk' M=N±l E' kk' 
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where n,m,k, and k' are indices of the eigenstates of the DWCNT QD Hamiltonian. The Redfield tensors have the 
form 



/ ALE'J 

-Ei/ _ -p{+)E'M En I -p(-)-E^j Bjv frn\ 
^nmkk' ~Z_^'^lk'mnk '^'^ I k'mnk ' 
/ 

and the matrix elements of the electron operators are 

i^U^))Z """^^ = (N, k |*L(x)| N + 1, m) 

with the states |N, k) and |N + 1, m) having energy En, E'j^^-^ and particle number N, N + 1, respectively. Such 
matrix elements are calculated in analytic form in App. [A) The transition rates depend on the properties of the 
contacts between the leads and the DWCNT QD^^. Here, we assume that the contacts are very simple. They do not 
mix the electrons in the different branches and the couplings between the leads and the DWCNT do not depend on 
either the wave vectors or the spins of the tunneling electrons. Then the transition rates depend on the energy of 
the tunneling electrons because of the matrix elements of the electron operators, Eq. (jA3p and have the forms (the 
derivation of these expressions is given in App. IB|) . 

riVm-j"'' - E '-Y^d^^^f^'O E n n n n 

a faF q>Qq'>0 jSli j'S'C' (53) 

X FiX^^^fg{ui),k'jg^q,mjSiq)F*{Xi~J/{ui),nj,s'eq'^kj,^ 

where the constants 7;^ describes the coupling strengths between the bonding/anti-bonding state a in c-DWCNTs 
or between the shell a in i-DWCNTs and the leads / and ui = 0, L ior I = s,d. The vector Bafa- denotes a state with 
one particle in the branch ara. The function F{\,ra,m') is defined in Eq. (jA4p and the parameters A;j^^^^(u)'s are 
defined in Eq. (|A5|) . The eigenstates involved are 

|fc')=|N,k'), |m)=|N + l,m), 

|n)=|N+l,n), |fc)=|N,k). 



Similarly, the expressions for the remaining tunneling rates are 

lie 

h 



rit'i"""-^ = E if 5.(^0(1 " /(^o) E ^N-e..„,N-i n n n n 



faF q>Oq'>0 jSi j'S'C (54) 

X ^Siq.'mjs^g)F{X^JJ/{ui),nj,s'^'g',kj,s'^'q') 



with the eigenstates 



|fc')=|N,k'), |m)-|N-l,m), 
|n)=|N-l,n), |fc)=|N,k). 

In the linear transport regime, only the following tunneling rates between the ground states with N and 7V±1 electrons 
are needed, which have very simple expressions, 

rltlfr""^ =J2Y.^9i{ei)m) (55) 
with the eigenstates |A;) =|N,0) and |n) =|N+ 1,0), and 

r|Zr"-^ =EEx^'(^')(l-/^^'))'5N-e„...N-l (56) 
a faF 



with the eigenstates |fc) =|N, 0) and \n) =|N — 1,0). We are only interested in the properties of the system in the 
stationary state, which can be obtained by solving the Eq. ([501) with the left hand side set to be zero. 
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B. Calculation of the current 



The current can be calculated by using the tunneling rates between the DWCNT QD and the leads. The current 
measured in experiments is the current in one lead, which can be calculated as 

/; = 62^(6; -Gj ), (57) 

TV 

where Q^^^^^ arc the tunneling rates between the QD and the lead I when the particle number in the DWCNT QD 
changes from to ± 1. The tunneling rates are related to the transition rates and the reduced density matrix as 

r^N^N±l /'-p(+)£^iv -E^^i i^En , /,-Eiv-p(-)-^« £n±i\ /cq\ 

^ 2^ l,nnk Pkn + Pnk ^ l,kjj7i ) ■ l&^i 

E,E' nkj 

After substituting Eq. ((58)) into Eq. ([57)) . the current can be expressed in terms of the transition rates and the elements 
of the reduced density matrix as 

— ^ \^ l,njjk ^ l,njjk J Pkn + Lkjjn ^ Lkjjn J Pnk ■ ) 

N.E.E' 



IV. LINEAR AND NONLINEAR TRANSPORT 



After having obtained the energy spectrum and the eigenstates of the DWCNT QD system, we can calculate the 
transition rates, Eqs. ([53]) and (|54p and use the Bloch-Redfield equations for the reduced density matrix to calculate the 
transport properties of the system. Here we present the calculated results of both linear and nonlinear conductances. 



A. Linear conductance 



In the linear transport regime, i.e., |eVb| ^ fc^T ^ Eq, where VJ, is the applied bias, only the ground states with N 
and A^ + 1 electrons are involved in the transport. In this case, the equations for the diagonal elements and the off- 
diagonal elements of the reduced density matrix are decoupled from each other and we only have to take into account 
the diagonal elements of the ground states with a certain electron number, which are the occupation probabilities. 
The stationary occupation probability of the ground state with A^ electrons is given as^^ 

p _ Y.la 7;a(l - /(£i))C^+l.Ar . 

Eio7/c/(e/)C^,^+i+7;a(l-/(£0)C^+i,Tv' 

where jy+i the number of permitted ground states with A^ + 1 particles when one electron is added to a ground 
state with A^ particles and this electron is added to the bonding/anti-bonding state a in c-DWCNTs or to the shell 
a in i-DWCNTs. We define the energy £i = eVi — AE and the energy difference AE — E^ — E^^-^ — jig, where fjLg is 
the electrochemical potential in the gate. The linear conductance is then given as 

h El. lU{-AE)C%,j,^, + 7,.(1 - /(-Ai?))C^^,_^ (1 + e-P^^Y ' ^ 

where we assume that the bias is symmetrically applied to the source and drain leads, that is, —Vs — Vd — Vb/2. The 
maximum value of the linear conductance is 

_ 2e^/3 X]q 7sQ7rfaC'^,Ar+iC'w'+i,JV 

Lrmax — — 1 (OZj 

Ula 7/Q(C^,Ar+l + '^N + I.n) + Ha'^N.N+l) iUla 7iQC'w+l,Ar) 

and the maxima of the conductance as a function of are —^ig — -Ejy+i ~ + ^-E'max, wherei^^ii^i^i^ 

AE^,^ = ^ In V ^ • 
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The conductance peak occurs whenever an electron is added or removed from the DWCNT QD by changing the 
electrochemical potential in the gate. At zero temperature, from Eq. (|63p ASmax vanishes and the conductance peak 
occurs when the electrochemical potential of the gate satisfies the following condition, 

-fig = E%_^_^ - E% = flN. 

Therefore, at zero temperature the addition energy S/in is given by 

SfiN = Imw - = \E%+, + E%_, - 2E%\. 

For a c-DWCNT QD system, electrons can tunnel into both shells because of nonzero intershell couplings. Hence 
there is an 8-electron periodicity of the conductance peak distances, which are 

Sfii = = Sfi5 = 6fj.7 ^W^^+, (64) 

6fi2 = = 2 min( A, C)eo + W++ , (65) 

(5/i4 = 2|A - Cko + W++, (66) 

,5^8 -eo-2(A + C)eo + Woy. (67) 

Here, we use the relation, W(^'^ — Wgo^ = Wqq^ /2, in c-DWCNTs (cf. Eq. ((26l) ). On the other hand, electrons can 
only tunnel into the outer shell in an i-DWCNT QD system because the contacts are deposited onto the outer shell 
and the intershell couplings vanish. Therefore, there is a 4-electron periodicity of the conductance peak distance like 
in a SWCNT QD system, which arc 

Sfii = dfi3 = W++, (68) 
Sfi2=2A+eo + W++, (69) 
dfii^eo-2A+eo + W++. (70) 

Because electrons tunnel only into the outer shell with the interaction strength WqJ,^, the addition energy S/in does 
not depend either on the interaction strength in the inner shell W^oo~ nor on the intrashell interaction strength Wq{,^ . 
The calculated linear conductances of DWCNT QDs of different configurations are shown in Figs. [5] and [HI In a 
c-DWCNT QD, the intraband interaction strengths W^'^ and Wqq~ are the same while the interband interaction 
strength W^qo^ twice as large (cf. Eq. ((26|) ) . However, in an i-DWCNT the interaction strength in the inner shell 
Wqq~ is the strongest because of the smaller inner shell radius and one has Wqq~ > W^'' > W^'^ . The shapes of the 
conductance peaks strongly depend on the mismatch parameters, i.e., A and ( in c-DWCNTs and A± in i-DWCNTs. 
For zero mismatch parameters, the quantities and C^+i n+i ~ 3, 2, 1 and C'^_^_l jy = 1, 2, 3, 4 for 

TV" = Am, Am + I, Am -\- 2, Am -\- 3 with an integer m, where N°' is the electron number either in the bonding/anti- 
bonding state a in c-DWCNTs or in the shell a in i-DWCNTs. Therefore, according to Eq. (p^ one can find 
that the conductance peak heights show an 8-electron periodicity in a c-DWCNT QD (cf. Fig. [5^a)) because both 
bonding and anti-bonding states contribute to the electron transport. However, there is a 4-electron periodicity in an 
i-DWCNT QD (cf. Fig. [SI|a)) because only the outer shell contributes. If the mismatch parameters are nonzero, we 
find C^_^+i = 2, 1, 2, 1 and C^+i_^ = 1, 2, 1, 2 for iV" = 4m, Am +l,Am + 2, Am + 3. Therefore, all the conductance 
peaks have the same heights (cf. Fig.El^b) andEUb)). However, the distance between two conductance peaks, i.e., the 
addition energy, always shows an 8-electron periodicity in c-DWCNT QDs and the 4-electron periodicity in i-DWCNTs 
as shown in the Figs. [5] and [6l 



B. Nonlinear conductances 



When higher bias is applied, i.e., |eVb| > £o ^ ^bT, we can only solve the Bloch-Redfield equations numerically. 
For elastic tunneling process, we have to include the coherences between the states with same particle number N but 
with different bosonic excitations m^^'^^. Because of the large number of degenerate bosonic excitations, the rank 
of the reduced density matrix increases very fast as the applied bias increases, which causes a very long computing 
time to solve the equations. On the other hand, these coherences can be ignored in an inelastic tunneling process, in 
which the QD system will be restored to the equilibrium states before the next tunneling process. Only the diagonal 
elements in the reduced density matrix are nonzero and they obey the Boltzmann distribution as 
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FIG. 5: (Color online) Calculated linear conductances as a function of the gate electrochemical potential in commensurate 
double- walled carbon nanotube (c-DWCNT) quantum dot (QD) systems with different parameters, (a): A = ( = 0.0, 
Wix)"'' ~ Wf^~ — Wq^~ /2 — S.Oeo, and fcsT — 0.025£o, where the level spacing eo is used as the unit of energy. The coupling 
strengths are 7s± — 'yd± = O.Oleo. (b): A = 0.2 and C, = 0.3. The remaining parameters are the same as those in (a). In both 
cases, the linear conductances in c-DWCNTs show an 8-electron periodicity. In the case of zero mismatch parameters shown in 
(a), an 8-electron periodicity of the heights of the conductance peaks also occurs. For finite mismatch shown in (b), the peak 
heights are equal but the 8-electron periodicity of the addition energies, i.e., the peak distances, remains, (to emphasize this 
we assign to each ground of eight peaks different colors). 
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FIG. 6: (Color online) Calculated linear conductances as a function of the gate electrochemical potential in incommensurate 
double-walled carbon nanotube (i-DWCNT) quantum dot (QD) systems with different parameters, (a) A+ = A_ = 0.0, 
Wqq^ = 5eo, M^oo~ ~ 6.0eo, W^oo~ ~ 5-5eo, and ksT = 0.025eo, where the level spacing eo is used as the unit of energy. The 
coupling strengths are 'ys+ ~ 7d+ ~ O.Oleo and 7s- = 7d- = 0. (b): A-i- = 0.2 and A_ — 0.3. The remaining parameters are 
the same as those in (a). The linear conductances in i-DWCNTs show a 4-electron periodicity. In the absence of mismatch 
shown in (a), the 4-electron periodicity is observed also in the peak heights while it is no longer observed at finite mismatch 
shown in (b). However, the addition energy, i.e., the peak distance, shows a 4-electron periodicity in both cases, (to emphasize 
this we assign to each ground of eight peaks different colors). 



where n and k are indices of the eigenstates of the DWCNT QD Hamiltonian and VN{t) is the probability of finding 
N electrons in the QD. Instead of solving the Bloch-Redfield equations directly, we can solve the equation of motion 
for the probability VN{t), 

jt'PN{t) = - E 0^^''+ E 0''^"^' (71) 

l,M=N±l l,M=N±l 

where the tunneling rate is defined in Eq. (j58p and can now be expressed in terms of VN{i) as 

^VN{t) 2^ , \X.^l,n33n +^l,njjn )■ (72) 

E,E' l^k ^ nj 

The number of the equations reduces significantly and the equations can be solved quite fast. In Fig. [71 we show 
the calculated stability diagram of a DWCNT QD system in an inelastic tunneling process. The size of the Coulomb 
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FIG. 7: (Color online) Calculated stability diagrams of double-walled carbon nanotube (DWCNT) quantum dot (QD) systems, 
(a) Stability diagram of a commensurate DWCNT(c-DWCNT) QD. The parameters are A = 0.2, ( = 0.3, W^o+ = Woq" = 
Wq^~ /2 = S.Oeo and fcsT = 0.05eo, where the level spacing eo is used as the unit of energy. We use 7a± = 7tj± = O.Oleo for the 
coupling strengths between the leads and the DWCNT QD. (b) Stability diagram of an incommensurate DWCNT(i-DWCNT) 
QD. The parameters are A+ = 0.2, A_ = 0.3, W^^'^ = 5eo, Wqq~ = 6.0eo, W^oo~ = 5-5eo and ksT = O.OSeo. The coupling 
strengths are 7^+ = 7^-1- = O.Oleo and 73- = yd- = 0. The size of the Coulomb diamonds shows an 8-electron periodicity in 
c-DWCNT QDs as shown in (a) while it shows a 4-electron periodicity in i-DWCNT QDs as shown in (b). 

diamonds shows an 8-electron periodicity in c-DWCNT QDs and an 4-electron periodicity in i-DWCNT QDs. The 
excitation lines are also shown in Fig. [71 which contain contributions of both fermionic (cf. Eqs. (PT]) and 1321)) and 
bosonic excitations (cf. Eq. (22)). There are more excitation lines in c-DWCNT QDs than in i-DWCNT QDs because 
of the larger number of the ground states of c-DWCNTs. The stabiUty diagram of an i-DWCNT QD looks quite similar 
to that of a SWCNT QD, which shows also a 4-electron periodicity. However, the configuration of the excitation lines 
of the two cases are different, because the excitation spectrum in i-DWCNTs contains an extra contribution from the 
Coulomb interaction due to the electrons in the inner shell. 



V. CONCLUSIONS 

In this paper, we derive the energy spectrum of both finite length c-DWCNTs and i-DWCNTs with open boundary 
conditions. Then we calculate the transport properties of the DWCNT QD system by solving the Bloch-Redfield 
equations for the reduced density matrix of the QD systems. Because the contacts are usually deposited on the 
outer shell and the intershell coupling depends on the chiralities of the two shells, we find an 8-electron periodicity 
of the linear conductance peak distances in c-DWCNTs but a 4-electron periodicity in i-DWCNTs. The peak heights 
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strongly depend on the degeneracies of the ground states. By including both fermionic and bosonic excitations, we 
also calculate the stability diagrams of QD systems with both c-DWCNTs and i-DWCNTs in an inelastic tunneling 
process. The periodicity of the Coulomb diamond sizes depends on the number of the shells contributing to the 
electron transport. Therefore, the 4-electron periodicity in a MWCNT QD measured in the experiments in Ref. '2l| 
may be due to the fact that only the outermost metallic shell was involved in the electron transport. Similarly, the 
4-electron periodicity in i-DWCNTs is because of the negligible intcrshcU coupling and there will be an 8-electron 
periodicity if a large intershell coupling is caused for example by contacts. Therefore, it's necessary to use properly 
prepared contacts in order to observe the different periodicities of Coulomb blockade oscillations in different types of 
DWCNTs in the experiments. 
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APPENDIX A: MATRIX ELEMENTS OF THE ELECTRON OPERATORS 



In this appendix, we calculate the matrix elements of the electron operators Eqs. (1^^ and (P7)) in the basis of 
the eigenstates, Eq. of the Hamiltonian Hqd and the results are used in Sec. IIIII The matrix elements of the 

electron operator are calculated by using the relations 



(N,m|*„a(r)|N',m') = VI^ sgn(F)(p„,gn(_F)i._F(r)(N, m|V'„f<TF(w)|N', m'). 



rF 



(Al) 



The ID electron operator i/jafFa can be expressed in terms of the bosonic operators introduced above as^i^^ 

i;^,p^{u) = !?£^f5£2f£4^e'*--(")+^*— (A2) 
VI - e-'^'^/^ 

where a is an infinitesimal positive number used to avoid the divergence in the long wave length limit and the operator 
rj^f^ is the Klein factor, which destroys a particle in the branch ara when acting on the eigenstates of the DWCNT 
Hamiltonian, 

r;„,JN,m> = (-l)^?^'r^ ^.|N - e„,., m), 

where we use the convention i — +R t, +L t, +R i, +L |, —R f , —L t, — i? J,, —L | = 1, 2, 3, 4, 5, 6, 7, 8 and the vector 
Gafo- denotes a state where there is only one particle in the branch afa. The notation X]r=i ^ means that the sum 
runs over all the state from 1 to i — afa — 1 with ara fixed by the unit vector e^fa- The phase factor is 

K -r. (lA - J^sg,n(F)(ssrn(f)No,f^ + A^)u 

^^arF(7\^) — fKT ' 
V 2L 



where A± = A for the c-DWCNTs. The field operator <j>l^fp„ is given as 



^ p-aq/2 

irh — > - p*sgn(r_F)t3u. 

^VafFa-y^) — / J . ^ 'Jasgn{r)qa 

q>0 V 



Therefore, the matrix elements of the ID electron operator have the forroi^ 

Vl - e-"'^/ 



n 40 



where the function F can be expressed in terms of the Laguerre polynomials L 

F{X,m,m') = ^^^^L™— ~"™(|An(e(TO'-m)A'"'-™ + e(m-mO(-A*)'"-™^ (A4) 

with TOmax = max(TO, m') and m^i^ = min(m, m'). 0(x) is the Heaviside step function and the parameters A's are 
given by 

\jS^q 
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APPENDIX B: EXPRESSIONS OF THE TUNNELING RATES 



In this appendix, we give a derivation of the expressions of the tunnehng rates, Eqs. (|53|) and (|54p . In general, the 
expressions of the tunnehng rates F's are given bj^, 



+ 1 
nk 



r!JS/"' ^x.y(*„,.,) ";-(*L(.)): 

X r dt' Ti^,{^,y,t')e^^-i^'«^-~^-y , (Bl) 
Jo 

X r dt'£iU^,y,t')e^iiE'.-.-E.)t' ^ (B2) 
Jo 

where t' — t — ii. The two functions are defined as 

(x, y, t') = Ti^ (x)T,; (y) (x)<f>L (y, -t'))^^ 

q 

^,„,(x,y,t') =T;Jx)T,„( «ey) ($t^(x)$,^(y, -t'))^^^ 

= T;jx)T,„(y) / de5,(e)/(£)E</.;(x)^q(y)e^(^-^^')*', (B4) 

q 

where (/((e) is the density of states in lead I, Vi is the voltage in the lead I, and /(e) is the Fermi distribution function. 
In the following we shall derive the expressions for the tunneling rates rjj,/^^^ for an i-DWCNT QD system as 
an example and the expressions for the other tunneling rates can be obtained by the same method. By substituting 
Eqs. dm and ^ into Eq. (|Bl]), we have 

rife^^fe''""^ = xEEE«g^^(^^') //rfxdy ^„,g„(^),^(x)^:,g„(^,),,^,(y) 

Qcr FF' 

X gi{e,)j{ei)Tt^{^)Ti^{y) E 'U^)Kky) (B5) 



where e; = eVi — -Eat — E'at+i and m; = 0, L for I = s,d. We ignore the slow oscillations of the ID electron operators 

IpafcrF (""0 ) ( raf'crF' ) ) 

independent of the positions. The integrals over e and t' are carried out by using 



J degiie) dt'e^^^'"^^'' = ^hg{E) ± iKP J j^de 



with V denotes the Cauchy principal value. We assume that the width of the lead energy band is infinite and that the 
lead density of states giie) is constant. Hence the Cauchy principal value is zero. Let's focus on the part depending 
on the position in Eq. (IBSp . namely, 

/ = / / dxdy Vasg„(F)fF(x)¥'Lg„(F')?'f '(y)^ia(x)r/a(y) E "^^q W'/'q (y ) • (^6) 

q 

Because the Bloch waves </3Qsgn(F)fF from Eq. ^ are largely localized around the carbon atoms and on the length 
scale of the Pz orbitals all the other quantities in Eq. (|B6|) are slowly varying, we can rewrite the two integrals as two 
sums over the positions of the carbon atoms and the Eq. (jB6p becomes 



^ = ]^ E E e^'''''e'^'''''/apsg„(F)fF/:j,.gn(F')f'F'7^/;(xR,p)Tz„(yR.,,0 E '/';(xR.p)<^q(yR'.P'), (B7) 

" Rp R'p' q 
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where the constant C denotes the integration over the pz orbitals. Because the leads are described by 3D Fermi gases, 
the wave functions (/i;q(x) are the plane waves, 



with the volume of the gas Vi. The sum over q can be performed 

47r sin(|q||x-y| 



^0;(x)0q(y) 



|q||x-y| 



which is peaked around x = y. Because of the large Fermi energy in the leads, the above expression can be approxi- 
mated by two Kronecker J's, 

q 

Therefore, Eq. (|B7[) becomes 

I irr X! "^''^"^ '^''^Iapsgn{F)fFfapst!,n(F')f'F'Tla{^Il,p)Tla{y-R.,p)- (B8) 
Rp 

Because of the fast oscillating phase e*^*^"^ the quantity / is nonzero only if F = i^' and in turn the sum over p 
can be carried out as 



/apsgn(_F)rF/Qpsgn(F)f'F — ^ff', 

P 

which can be easily verified by using the explicit expressions Eqs ([7]) and ([8]) . Eq. (|B8|1 becomes 

^ - ^ E Til{^B^,p)Tio.{yii,j,)SFF,6,f,. (B9) 

Rp 



By substituting Eq. (|B9[) and Eq. (|A3|) of the matrix elements of the electron operators into Eq. (jB5|l . we finally 
obtain the expression of the tunneling rates, 

rit'i" ""^^ = E ^fa^i^om E ^N+e„..,N+i n n n n 

a faF q>Oq'>0 jSU'S'e (BIO) 

X ^(^afi^("/). kjSiq^rnjSiq)F* (Xl^-J^^ (Ul), n^'S'eq'^kyS'i'q')^ 

where the vector eaf^ denotes a state with one particle in the branch ara. The function F{X, m, m') is given in 
Eq. (jA4p and the parameters X-^^^f's are defined in Eq. (jASp . The constants jia describes the coupling strengths 
between the shell a in i-DWCNTs and the leads / having the form 

|T';a(xR. 



The four eigenstates are 



Rp 



|fc')=|N,k'), |m)=|N + l,m), 

|n)=|N + l,n), |fc)=|N,k). 

Similarly, the expressions for the other tunneling rates are 

r[v^:,^-' - E ^9i{ei){i - /(so) E '5N-e„.„, N-i n n n n 

a faF q>Oq'>0 jSli j'S'e (Bll) 



X P*iK,iqM^ kMq^^35iq)F{Kfaq M , ^j' S' ^' q' , kj, S' ^ q') 
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with the eigenstates 

|fc')-|N,k'), |m) =|N-l,m), 

|n)=|N-l,n), |fc)-|N,k). 

From the expressions of the tunnehng rate, Eqs. (|B10|) and (jBll|) . we can see clearly that the contacts do not mix 
the electrons in different branches. 
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